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Abstract 

Based on a modified expression of the rate of the convective constraint release, we present a 
new contour-variable model of constitutive equations in which the non-uniform segmental stretch 
and the non-Gaussian chain statistical treatment of the single chain are considered to describe the 
polymer chain dynamics and the rheological behavior of an entangled system composed of linear 
polymer chains. The constitutive equations are solved numerically in the cases of steady shear and 
transient start-up of steady shear. The results indicate that the orientation and stretch, as well as 
the tube survival probability, have strong dependence on the chain contour variable, especially in 
the high-shear-rate region. However, the inclusion of the non-uniform features in the constitutive 
models has little modification comparing with the uniform models in determining the rheological 
properties both qualitatively and quantitatively. 
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I. INTRODUCTION 



Constitutive equations in polymer melts or concentrated solutions are mathematical rela- 
tionships between the stresses and the external flow conditions. They represent the inherent 
properties of the polymer system, and should be derived from the knowledge of polymer 
chain structures, configurations, interactions and polymer dynamics. In practice, a consti- 
tutive model is needed to idealize these microscopic characteristics of polymer chains, due to 
the complex interactions of this many-chain system (Larson, 1988). As polymer melts own 
some characteristics of universality, that is to say, some physical properties of polymer melts 
are independent on the microscopic chemical structures, we can coarse grain a real polymer 
chain as a smooth thread connected by segments (Watanabe, 1999). This coarse-grained 
chain, just like the real chain, has an enormous number of configurations and can be usefully 
described statistically as Gaussian chain or non-Gaussian chain. In fact, the dynamics in our 
study is the segmental dynamics, which is based on a spatial scale of a chain segment and 
those chains are only distinguished by their lengths and abilities of extension. The interac- 
tions between chains are intricate. However, the attractive interactions between segments 
are tended to screen out by the excluded volume effect in polymer melts or concentrated 
solutions (de Gennes, 1979), leaving the topological constraints to play a dominant role in 
such a system at equilibrium. A basic understanding of chain dynamics is the key-point 
to develop constitutive equations for melts or concentrated solutions of flexible polymers. 
The proposition of the idea of 'reptative' motion by de Gennes (1971) is a milestone in 
the development of polymer dynamics of the concentrated system. In fact, in a flow the 
long linear polymer chains become oriented and stretched, which generates internal stress 
as the response to external disturbance. However, the orientation and stretch of chains will 
relax in the course of time due to the motion of non-crosslinked polymer chains. The major 
mechanisms of relaxation concerned are: reptation (Doi and Edwards, 1978a, b, c, 1979), 
convection of segments along primitive chains (Doi, 1980; Marrucci, 1986), fluctuations of 
the contour length (Doi and Edwards, 1986; Pearson and Helfand, 1984; Mead et al., 1998), 
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reptative constraint release (Marrucci, 1985; des Cloizeaux, 1988, 1990; Leygue et al., 2005, 
2006a, b), and convective constraint release (CCR) (Marrucci, 1996; Ianniruberto and Mar- 
rucci, 1996; Milner et al., 2001; Graham et al., 2003). The first three kinds of relaxation 
are caused by the motion of the primitive chains, which can happen only at the ends; while 
the other two are induced by the motion of surrounding chains, which can occur not only 
at their ends but also in some other place along the primitive chains. 

The optical-stress law guarantees a linear relationship between the stress and the orien- 
tational anisotropy of segments for a highly crosslinked rubber system (Treloar, 1975). The 
subsequent studies also show its validity for non-crosslinked polymer system (Janeschitz- 
Kriegl, 1975). The segments in the crosslinked rubber refer to chain portions between two 
permanently tethered crosslink points. The orientations of these segments are uniform. The 
segments in the non-crosslinked polymer system refer to entanglement segments. The ori- 
entations of these segments are non-uniform, as the segments can move freely. Under fast 
flow the polymer chains not only become oriented but also significantly stretched. The con- 
tribution from segmental stretch to stress is proportional to the square of stretch ratio for 
Gaussian chains (Doi and Edwards, 1986), while this relationship is much more complex 
for non-Gaussian chains. In the original Doi and Edwards (DE) model, due to the faster 
relaxation of segmental stretch, the chain is assumed to remain its equilibrium length all the 
time. The stress contribution comes only from segmental orientation. DE theory predicts 
a maximum in shear stress a xy with respect to the shear rate, and followed by a sharp de- 
crease asymptotically as 7 -0 ' 5 , which has not been observed by experiments. The neglect of 
segmental stretch results in a monotonic increase in the first normal stress difference in the 
start-up of a simple shear flow, however, an overshoot appears in experiment. In order to 
improve this deficiency of the original DE model, several 'chain stretch' models are proposed. 
In the uniform stretch model, it is assumed a uniform segmental stretch along the chain, 
such as Larson's 'partially convected strand' model (1984) and also the one proposed by 
Pearson et al., (1989). The non-uniform segmental stretch model is proposed by Marrucci 
and Grizzuti (1988). The segmental stretch models are called as DEMG model (Mead and 
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Leal, 1995a, b). The extension of chain segments creates some segment (or saying pseudo 
defect (de Gennes, 1971)) sources along the chain, which is equivalent to reduce the diffu- 
sion coefficient. The inclusion of chain stretch also changes the feature in the start-up of 
shearing. A serious deficiency in all these models is that in the terminal region excessive 
shear thinning is inevitable due to the lacking of an extra relaxation mechanism to increase 
the rate of creating new segments, which is equal to the rate of annihilating old segments, 
or equivalently saying to avoid high orientation of the chain to the direction of flow. The 
Mead-Larson-Doi (MLD) model incorporating CCR to the DEMG model qualitatively im- 
proves the description of the phenomenon of excessive shear thinning and other rheological 
properties in steady state and transient shearing flow (Mead et al., 1998). Recently Pat- 
tamaprom and Larson (2001) proposed a toy MLD model by making an extension of the 
MLD model in a simple way and adding reptative constraint release to MLD equations. In 
all these models, the rate of CCR is evaluated at the tube end, which is a uniform rate of 
CCR. 

CCR represents the ability of the flow to convect or release entanglements, and is deter- 
mined by segmental orientation and stretch. As the segmental orientation and stretch in the 
non-crosslinked polymer melts are non-uniform, the CCR is non-uniform. The rate of CCR 
is a function of both the contour variable and time. We assume that when the entanglement 
is convected away from the primitive chain segment to some other place, the previously 
formed environment has changed and this entanglement then becomes released. Based on 
this idea, we develop a new contour- variable model including non-uniform segmental stretch, 
non-uniform rate of CCR and non-Gaussian chain to describe the rheological behavior of 
linear chain entangled system in simple shearing and transient flows. However, we dropped 
the reptative constraint release term in the present work as we focus on the high-shear-rate 
region, where the segmental orientation and stretch are obvious. In the low-shear-rate re- 
gion, the reptative constraint release occurs on a time scale of the reptation time of the whole 
chain and one event only relaxes a small part of the chain (Graham et al., 2003); while in 
the high-shear-rate region, the CCR dominates. Thus in both of these two cases it is valid to 
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drop the reptative constraint release term. While in some middle region, some quantitative 
deviations occur if we drop the reptative constraint release term. We also neglected the 
contribution from the fluctuations of the contour variable, as the problem is considered in 
the mean-field level. 

This paper is organized as follows. In section II we develop a contour- variable dependent 
rate of CCR to describe the flow induced constraint release and construct the new equations 
which the segmental stretch and the tube survival probability satisfy, respectively. We 
promote our contour-variable model of constitutive equations in the end of this section. In 
section III we give the results and discussion, in which we focus on the simple steady shear 
flow and the start-up of steady shear flow. In section IV, we give the conclusions. The 
process and methods used for the calculation of constitutive equations are shown in the 
Appendix. 

II. THE NEW CONTOUR- VARIABLE MODEL 
A. Rate of CCR 

CCR represents the entanglements relaxation mechanics caused by the convective motion 
of surrounding chains. It is not important when the concentrated system is at rest or 
undergoes a slow flow. However, when the velocity of the flow is comparable to the inverse 
of reptation time of the primitive chain, CCR will play a key role in determining the dynamics 
of the system. Taking into account the CCR mechanics, Ianniruberto and Marrucci (1996) 
proposed a contour location independent rate of CCR. In this model, the rate of CCR is only 
caused by the retraction of the primitive chain ends. So the rate of CCR is proportional 
to the velocity gradient of the movement of the primitive chain ends. However, not all 
convections or deformations will release constraint. In order to exclude the case of affine 
deformations, in which no constraint release occurs, Mead et al. gives the rate of CCR as 
the difference between the rate of convection of the entanglements and the rate of primitive 
chain retraction (Mead et al., 1998; Viovy et al., 1983). They only evaluate the rate of CCR 
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at the chain ends, since they think the topological constraint can be released only by the 
motion of chain ends (Ianniruberto and Marrucci, 1996; Marrucci et al., 2001). (We call it 
the end-relaxation mechanism in the follows.) 

In the present work, we focus on the non-uniform behavior of CCR. Let us consider a 
primitive chain segment confined by a tube segment which was created at past time t', as 
called the t'-segment. The deformation of this chain segment will be released when the en- 
tanglements forming the tube segments are convected away by the flow. It is an irreversible 
process. This i'-segment will never be created again. When this happens, some hidden en- 
tanglements become active (Ianniruberto and Marrucci, 2000), which instantaneously makes 
the released part (e.g. the t'-segment) become a part of a new tube. Thus we can say that 
the rate of creating a new segment is the same as that of annihilating an old segment. The 
contribution to the stress for a specified tube only comes from the part that remain unre- 
leased. The released part of the original tube, which becomes a part of a new tube, will 
still have contributions to the stress. However, the stress should be calculated from the 
new tube. Define Sq as the equilibrium contour variable, which runs from —L /2 to L /2, 
where L is the equilibrium length of the primitive chain. When the chain is stretched, its 
primitive length L can be larger than L . We then define a local stretch function s(so,t) 
which takes the value from —L/2 to L/2, to describe the segmental stretch. We can define 
q(so, t) = ds/dso as a local strain. It describes the local extension. The rate of CCR k(so, t) 
can be written as the following form 

k(s , t) = K: S(s , t)q(s , t) - 1 t_ (-q 

q{s ,t) dt 

where n is the velocity gradient tensor of the flow, which, we suppose, is homogeneous, 
although in general it is not; S(so,£) is the orientational tensor, S(s ,t) = (u(sq, i)u(so, t)). 
Apparently, it is a local approximation form. The first term in Eq. (1) describes the rate of 
entanglements convection at contour variable sq and t. The second term is used to exclude 
the case of affine deformation, in which there are no environment changes with respect to the 
primitive chain and no CCR happens. As S(±Lq/2,£) = 0, q(±Lo/2,t) = 1, from Eq. (1) 
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we then find no CCR happens at the chain ends of the primitive chain. In fact the ends 
of the tube change instantaneously with the ends of the primitive chain during movement, 
there is no relative motion between the ends of the tube and the ends of the primitive chain. 
The primitive chain and the tube are intrinsically synonymous, we treat them as different 
things only for convenience in describing the motion of the chain considered and the motion 
of its surrounding chains. 

The difficulty of the flow-reversal problem has be pointed out (Wapperom and Keunings, 
2000; Ianniruberto and Marrucci, 2001). If the present model reduces to the uniform and 
single-relaxation-time version, this problem can be amended using the absolute value of 
|k : S| q in Eq. (1). In the present non-uniform case, we suppose that the above expression 
can amend this difficulty In the following we will not write the absolute value sign in the 
rate of CCR by considering the flow only in one direction. 



B. Segmental Stretch 

The dynamics of segmental stretch is described by the equation of motion, in which we 
neglect the acceleration term. It is derived from the force balance between the segmental 
extension force caused by the flow and the entropic elastic force. If the chain is considered 
to be non-Gaussian, this equation can be written as (Mead and Leal, 1995a, b) 

'd 2 s\ 



d , , . , 3pZD 
-s(s ,t) = {v{s ,t)} + 



dq 



(2) 



dt a 

Here (v(s Q ,t)) is the relative pre-averaged tangential velocity with respect to the center of 
chain 

(v(s ,t)) = k : / S(s' ,t)q(s ,t)ds' . (3) 

J 

The parameter Z = M/M e is the number of entanglements in a primitive chain, where M e 
is the molecular weight between two entanglements. Z = T^/3Tr (Doi and Edwards, 1986), 
where is the reptation time and Tr is the Rouse relaxation time. D defines the diffusion 
coefficient along the primitive path, and is related to the reptation time by = Lq/tt 2 D. 
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(3 is a coefficient related to non-Gaussian behavior, which is given by 

I -a 2 







a" 



(4) 



where a = Lo/L max , and L max is the maximum length the primitive chain can be stretched 
to. (3 ensures that the magnitude of a fictitious force acting along the primitive chain is 
equal to 3/c#T/a when there is no segmental extension (Doi and Edwards, 1986; Marrucci 
and Grizzuti, 1988). Here, a is the step length of the primitive chain and is a constant in 
the present model. A good discussion for its magnitude is given by Milner (2005). The 
inverse Langevin function C~ l (x) with the fractional extension x = aq is used to describe 
the tension of a non-Gaussian finitely extensible polymer segment (Treloar, 1975), which 
satisfies the following equation 



C{x) = coth(x) . 



x 



(5) 



Expanding £ 1 (x) in a Taylor series and using Pade approximation (Cohen, 1991), we can 
obtain 



C-\x) 



9 s 297 5 

= 3x + -x 6 H x b 

5 175 

3-x 2 
~ x 



X 



2 ' 



(6) 



We differentiate both sides of Eq. (2) with respect to so and obtain 



dq(s ,t) , 3ZD(3d£- 1 (aq)d 2 q 3ZD(3 
= k : S(s , t)q H ; tt^- + 



dt a dq 8sq 

Insert Eq. (6) into the above equation, we obtain 



a 



d 2 
dq 



2 C-\aq) 



dq s 
9s 0j 



(7) 



d( l c i ovn^ 3 + aV d 2 q 4a 3 q 3 + 12ag ( dq 

- = « : Sq + ^ + 3Z^a- 



;i-«V) 3 



dsr 



(8) 



Taking into account the effect of CCR and following the argument by Mead, Larson and 
Doi (1998), finally we obtain 



dq 
dt 



or7T ^ Q 3 + aV d 2 q QVn ^ 4a 3 g 3 + Ylaq ( dq 

3ZD P(T^)d4 + 3ZDf3a ( i- aV f 

+ .:S ? -i(.:S,-i|)( ? -l ) . 



(9) 



The parameter 1/2 comes from the fact that the magnitude of the rate of entanglements 
reduction due to CCR is two times as large as the rate of segmental retraction. The boundary 
condition and the initial condition are 

q(s ,t)\ So=±Lo/2 = 1 

and 

g(so,*)|t=o = 1, 

respectively. The boundary condition is specified by the untethered condition of the chain 
ends which can never be stretched. 

C. Probability of Tube Survival 

The process of the segmental orientation can be described by the tube survival probabil- 
ity G(s,t,t'), which means that a tube segment created at past time t' still can be found 
at location s and time t. In the original DE model, G(s, t, t') satisfies a diffusion equation. 
When the segmental stretch is taken into account, such as the DEMG model, a convective 
term has to be added to the DE model. This term accounts for the flow induced convection 
of the defects to the ends of the primitive chain and the elongation of the tube segments. 
However, the orientation relaxation induced by constraint release has not been considered. 
Without that it makes the prediction of the rheological properties by these models qualita- 
tively different from the experimental results. 

Based on our understanding for the effect of constraint release, we modify the CCR term 
in the equation of the probability of tube survival in MLD model and drop the term which 
comes from the reptative constraint release. Including the contribution from chain stretch, 
the equation for the probability of the tube survival is given by 

dG(s,t,t') _ ^ d 2 G(s,t,t>) dGjsXf) 

dt ds* {V[S,)) ds ' 1 j 

where s = s(so,t) is the tube strain function. The boundary condition and the initial 
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condition are 

G(s,t,t')\ s=±L/2 = 

and 

G(s,t,t% =t , = 1, 

respectively. It is convenient to transform the independent variables s and t to other inde- 
pendent variables so and t (Pearson et al., 1991). 

dG(s ,t,t') _ D d 2 G(s ,t,t') , ( Ddq 1, , lgg\ dG(s ,t,f) 

dt ~ q 2 ds 2 , + V g3 9sQ g W°'^ + gatJ dso ' 1 J 

where g = q(so,t) is the function of local strain. When we consider the CCR mechanism, a 
new term must be added to Eq. (11), which is 



ldq\ 



K:Sq---^\G(s ,t,t') 



Here, S is the orientational tensor. Then we obtain 



dG(s ,t,f) D d 2 G(s ,t,t') ( D dq 1 1 ds\ dG(s , t, if) 

-{v{s , t)) H — 



dt q 2 dsl \ q 3 ds Q q ' q dt J ds 

-'©H*-^)^''' ' (i2) 

where f(ds/dso) = 1/q is the switch function (Mead et al., 1998). When q is large, the 
contribution from the last term in the right hand side of Eq. (12) is much smaller. However, 
when q approaches to unity, this term plays an important role. This is arisen from the large 
difference in the time scale between the reptation time and the Rouse retraction time. It 
means that the relaxation of the tube orientation caused by constraint release of surrounding 
chains should start to happen just when the primitive chain nearly finish its retraction. The 
boundary condition is 

G(sq, t, t )|s =±L /2 = 0, 

and the initial condition is 

G(s Q ,t,t')\ t=t > = 1, 

Note that G(so,t,t') is the probability function, which has a range from to 1. 
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D. Constitutive Equations 



The stress is the sum of tensile force of the primitive chain projecting to some direction u 
by averaging over the conformations of primitive chains. Those primitive chains penetrate a 
unit area plane with its normal direction u. The tensile force acting on segment sq, along the 
non-Gaussian primitive chain at time t can be expressed by the inverse Langevin function 

FM = *jf C -r (13 ) 
When the Pade approximation is used, we have 



F(s ,t) = ^-P- ^, 14 

1 — a z q A 



where (3 has the same meaning as in Eq. (4), and is used to conform that when q — 1, 
F(s ,t) = "iksTja is equal to the fictitious force acting on every segment (Marrucci and 
Grizzuti, 1988). The stress can be written as 

c f L / 2 

a = Tt( dsF(so,t)u(s,t)u(s,t)) 
Jy J-l/2 

c r L / 2 

- T7 / ds(F(s , t))(u(s, i)u(s, t)) 

JM J-L/2 

15 G° N /-W2 o 3-aV 



(3- —q S(s ,t)ds , (15) 

4 L j-Lo/2 1 — a z q z 

where c is the number of polymers in a unit volume and N being the degree of polymerization. 
We decouple the average segmental tension from the average segmental orientation in the cal- 
culation of the first line to the second line. G° N is the plateau modulus, G° N = 4ckBTL /5Na. 
S(so,i) can be expressed by 

sko = J[jt , dG{s ^y\ m,t'))- (is) 

Here, Q is the DE strain tensor without the independent alignment approximation given by 

/ (E-u)(E-u) \ 1 
^ \ |E-u| / (|E.u|) ' (17) 

where E is the deformation tensor (Doi and Edwards, 1986). 
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III. RESULTS AND DISCUSSION 



The results reported in this paper will focus on steady and start-up of simple shear flow 

of monodisperse concentrated polymeric systems. We scale the equations for simplicity. The 

time variable is scaled by the reptation time, t/Td — > t, and the spatial variable is scaled by 

the length of primitive chain at equilibrium, Sq/L q — > sq. For convenience, we use the same 

symbol to denote the scaled quantities. Eqs. (9), (12) and (15) become 

dq 3Z(3 3 + aV d 2 q 3Z(3a 4a V + 12aq ( dq V 
dt ~ vr 2 {l-a 2 q 2 ) 2 d? + ~^ 2 (1-aY) 3 V^o/ 

+n:Sq- 1 -(^:S q -^y q -l), (18) 
dG(s ,t : t') 1 <9 2 G(s ,M') , ( I dq 1 lds\dG(s ,t,t') 

\V{Sq, t)) H 777 



dt n 2 q 2 8sq \ ir 2 q 3 dso q q dt J ds 



o 



^ 7^q 2 S(s ,t)ds , (20) 



■V* ? 3-a 2 q 2 ^ 
-i/2 1 — a 2 q 2 

respectively. Here, (v(so,t)) is given by Eq. (3), the relationship of which before scaling and 
after scaling is (v(sq, t))Td/ L — > (v(so, t)). The stress tensor has been scaled as 4a/15G° N — > 
a. S(sq, t) is given by Eq. (16). For the case of shear flow, k = ^e y e x . The shear viscosity is 
defined by r\ = <7 xy / A f, where a xy is the shear stress and the primary normal stress coefficient 
is given by \& = (<y xx Cyy^/if^i where u xx <Jyy is the first normal stress difference. The 
boundary conditions are 

g(s ,*)L=±i/2 = 1, 
G(so,t, OL=±i/2 = 0, 

and the initial conditions are 

q(s ,t)\t=o = 1, 
G(s Q ,t,t% =t , = 1, 

These equations are nonlinear and have to be solved numerically. Combining the finite dif- 
ferencing and the Newton-Raphson method, we can obtain the solutions of these equations. 
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The details are shown in the Appendix. The constitutive parameters in these equations 
include Z and a. Here, Z is the number of entanglements per chain and can be specified 
by fitting the steady state predictions of the model and the experimental data. A larger Z 
corresponds to a shorter Rouse relaxation time, a is the ratio between the equilibrium and 
the maximum length of the primitive chain, reflecting the extensibility of the chain. 

A. Steady state shear flow 

In this section, we show the dynamic properties of the chain segments and rheological 
behavior under steady shear flow. 

Fig. 1 gives the non-uniform segmental stretch of the primitive chain with different shear 
rates. The parameters are chosen as Z = 20 and a = which denote the Gaussian chain. 
Due to the untethered fact of the chain ends, the segments at chain ends remain unstretched. 
There exists a maximum value of stretch at the center. The extension of the segments 
equivalently reduces the diffusion coefficient by 1/q 2 in Eq. (19). So the equivalent reduction 
of diffusion is much more at the center. In addition to a small value of the rate of segmental 
convection around the center, we find that the rate of tube relaxation, or the rate of tube 
creation, is relatively small at the center. The cancellation of gradually increasing rate of 
convection and decreasing rate of effective diffusion from the center to ends induces a flat 
region of the rate of the tube relaxation around the center. According to the calculation 
from Eq. (16) under steady state shear, we find that around the center of the tube the 
value of S xy is relatively small. At the ends, where the segments of tube disappear at their 
creating time, S xy (±l/2,t) = 0. Therefore, there must be a maximum value of S xy in some 
location between the center and the ends. This can be seen from Fig. 2 when the shear rate 
is high enough. It is very interesting that when the shear rate is high enough, the maximum 
contribution to stress from segmental orientation comes neither from the ends of the chain 
nor from the center, but from some other location in between. Moreover, this location will 
move towards the ends when the shear rate grows higher. From the calculation we also 
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find that the segmental orientation is insensitive to the entanglement number of a primitive 
chain, denoted as Z which determines the segmental extension in steady state. Fig. 3 shows 
the differences of S xy between the present model and DEMG model for Z = 20 and Z = 50, 
respectively. The result shows that including of the CCR term will obviously enhance the 
value of S xy . Therefore, the magnitude of the stress in fast flow is seriously underestimated 
in the DEMG model. 

The comparisons of the shear stress and the steady state shear viscosity with different 
values of Z are given in Fig. 4. In the linear region, where the steady shear rate is smaller 
than the inverse of reptation time, both the shear stresses and steady state shear viscosity 
are independent of Z, because in these cases the shear rates are not high enough to stretch 
the chain segments. When the shear rate is high enough, the segments will be stretched to 
many times of their equilibrium length. The larger extension of the segment is, the more it 
contributes to the shear stress. On the other hand, the Rouse relaxation time Tr determines 
the magnitude of segmental extension. If Tr is larger, which corresponds to a smaller Z, the 
shear stress will be larger. Fig. 4 also shows the shear viscosity of the steady state which is 
defined as 77(7) = <y xy /^. When the shear rate is smaller than the inverse of reptation time, 
it behaves as Newtonian fluid. While in the high-shear-rate region it markedly depends on 
the shear rate. Different segmental stretch results in the differences of the first normal stress 
difference in the high-shear-rate region, which can be seen in Fig. 5. The same characteristics 
for the primary normal stress coefficients are also shown in Fig. 5. 

Up to now the chains are approximated as Gaussian chains {a = 0), which is valid if 
the stains are not too large. However, under fast flow the chains are largely stretched and 
approach to their limiting extensible values, the Gaussian treatment is no longer valid. In the 
following, we will use the non-Gaussian treatment of the single chain taking into account the 
finite extensibility of the chain (a ^ 0). Fig. 6 shows the comparisons of the shear stress and 
the first normal stress difference with respect to the shear rate with different extensibilities 
for a = 0, 0.2, 0.4 and 0.6, respectively. The parameter Z is chosen as Z = 20. The 
calculation also indicates that for a specified value of the shear rate segmental orientation is 
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not sensitive to segmental extensibility. From Fig. 6 it can be found that the curves of shear 
stress and first normal stress difference are independent of segmental extensibility in the 
small shear rate region where the stress contribution from segmental stretch can be ignored. 
In the high-shear-rate region, the differences of shear stress and the differences of the first 
normal stress caused by the differences of segmental extensibility are mainly contributed 
from segmental extension. The segmental orientation is only determined by the value of 
the shear rate, and is independent of the segmental stretch, because they are considered as 
independent processes. 

Fig. 6 also shows the comparisons between the predictions of the present model and the 
experimental data (Pattamaprom and Larson, 2001) of the stresses and the first normal 
stress differences under steady shear for different parameters of a; Fig. 7 shows the com- 
parisons with experimental data of the stresses and the first normal stress differences under 
steady shear for different parameters of Z. By fitting the experimental data, we find that 
incorporating the non-uniform features to the constitutive models does not affect Theolog- 
ical properties. Moreover, the non-uniform model gives us the detailed descriptions of the 
chain dynamics. In the present model there are two modifiable parameters of Z and a. By 
fitting the experimental data to the predictions of the present model, we can obtain the 
values of Z and a, respectively. The comparison between them will show you the relative 
abilities of the retraction and the extension of the chains. This means that, under the same 
magnitude of shear rate, different abilities of retraction or extension will result in different 
magnitudes of stresses, although the molecules composing the materials may have the same 
volume fractions and molecular weights. From these two figures, quantitative deviation from 
the experimental data exists, especially the first normal stress difference in the high-shear- 
rate region. These deviations may attribute to the "nonlocal" interactions between chains. 
By comparing with the uniform models by Pattamaprom and Larson (2001), we also find 
that the effect of non-uniform CCR on the stress and the normal stress is minor. The de- 
tails of segmental orientation and stretch are not important in determining the rheological 
properties in steady simple shear. 
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B. Start-up of Steady Shear Flow 

The transient behavior of chain stretch and orientation, as well as the rheological prop- 
erties, are discussed in this section. The steady shear is imposed to the system at t = 0. 

In Fig. 8 the segmental stretch varying with time at different shear rates are shown, 
and the constitutive parameters are chosen as Z = 20, a = 0.0. The contour variable is 
focused at s = 0. Before this segment reaches its steady value of extension, it passes a 
maximum value which is dependent on the shear rate, since other parameters are fixed. 
This maximum value will appear earlier when the shear rate increases. Fig. 9 gives the 
behavior of segmental orientation S xy (0,t) with the constitutive parameters chosen as the 
same in Fig. 8. We focus on the location sq = 0. If the shear rate is of the order of 1/Td 
or larger, the value of S xy (0,t) will also pass a maximum value before reaching its steady 
value. The appearance of the maximum value is not due to the fact of segmental stretch, 
but the existence of a maximum value in the strain Q xy with respect to the time. 

Fig. 10 and Fig. 11 show the transient behavior of the shear stress and the first normal 
stress difference, respectively, where overshoots are predicted in the present model. For the 
shear stress, when the shear rate is higher than the inverse of reptation time, the overshoot 
appears. The overshoot of the first normal stress difference appears at a much higher shear 
rate in the present model than that in other models. The reason is that the appearance of 
maximum in the first normal stress difference is determined by the evolvement of segmental 
stretch which should also have a maximum in the course of time. However, the rate of CCR 
we obtained contains a local strain q(s,t), which largely decreases the maximum value of 
segmental stretch and subsequently the stress. Thus the magnitude of the start-up shear 
rate, which can produce a overshoot, increases. 

Fig. 12 and Fig. 13 show the comparisons of the shear viscosities and the first normal 
stress differences with the experimental results, respectively. We choose the shear viscosi- 
ties instead of the shear stresses for convenience to compare with the experimental data 
(Pattamaprom and Larson, 2001), in which the material is sample one. The parameters in 
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these two figures are both chosen as Z = 20 and a = 0, 0.3, 0.6, respectively. The results 
indicate that only when the shear rate is high, i.e. 7 > 10/Tj, the differences caused by the 
different values of a are predicted by the present model. If the chains are less extensible, 
which corresponds to a large value of a, the segmental extension will be smaller under the 
same shear rate at the same time. This leads to smaller shear viscosities as well as the first 
normal stress differences. In fact, the segmental orientation is mainly determined by the 
magnitude of the start-up shear rate. When the shear rate is fixed the contribution from 
orientation to shear stress is specified. Thus, the magnitude of segmental stretch plays a 
critical role in determining the final value of the shear stress, and subsequently the shear 
viscosity. This is kept when the first normal stress difference is considered. By fitting the 
experimental data we find that the predictions of the present model have the same tendency 
with the experimental curves. On the other hand, there are large inconsistencies when the 
quantitative values are concerned, especially the magnitude of overshoots. Both the shear 
stress and the normal stress overshoots in the present model are very weak comparing to 
the experimental data. It is obvious that the chain stretch is not enough, which may arise 
from the inclusion of local approximation in CCR in the present model. 



IV. CONCLUSIONS 

In this paper, we develop a contour- variable model to describe the non-uniform behavior of 
segments under steady shear and start-up of steady shear, and the subsequent contributions 
to the shear stress and the first normal stress difference. The present model is a modification 
of the MLD model with a modified CCR term and non-Gaussian chain treatment of the 
single chain. We drop the contour length fluctuation contribution and neglect the reptative 
constraint release term, since both of them contribute less than the shear rate dependent 
CCR term in the high-shear-rate region. The numerical results predict strong dependencies 
of the segmental stretch and orientation on the contour variable. The maximum extension 
happens in the center of the chain due to the demand of symmetry. The fact that segmental 
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extension passes a maximum during the start-up of steady shear flow is the key reason 
why the first normal stress difference has a maximum. The segmental orientation also 
shows a maximum in the center when shear rate is small, and this is mainly attributed 
to the reptation. However, when the shear rate is high enough the segmental convection 
begins to work. A maximum of segmental orientation appears in some locations between 
the center and the end, where the stress contribution from orientation is the largest. The 
inclusion of CCR entirely promotes the magnitude of orientation of each segment rather than 
changes the distribution of the value of segmental orientation. The non-uniform segmental 
orientation and stretch subsequently contribute to rheological properties, which qualitatively 
agree with experimental data, although quantitative deviation still exists. When the shear 
rate is high, i.e. 7 ^> , the segmental strain is larger than 1. The difference of the 
segmental extensibilities is the main reason for the difference of rheological behavior. The 
same properties hold in the start-up of steady shear. This validates our pre- approximation 
of the independence of segmental stretch and orientation. 

A few approximations have been used in the present model. One is the local approx- 
imation of CCR. CCR happens once a tube segment is convected away from its original 
location, since the relative position to the primitive- chain segment confined by this tube 
segment changes. The convection of a segment may not be only determined by its local 
conditions, but all other segments elsewhere. In fact, the present model only predicts the 
overshoots in shear viscosity and first normal stress difference under the start-up of steady 
shear, but does not predict undershoots which have been observed in experiments (Huppler 
et al., 1967; Mhetar and Archer, 2000). It is well known that the appearance of over- 
shoots is due to segmental stretch. How about the undershoots? Pattamaprom and Larson 
(2001) predicted that the undershoots would appear if the dependence of the orientation 
and stretch on the contour variable was taken into account. The present model indicates 
that the inclusion of "local" non-uniform rate of CCR to the constitutive equations does 
not predict undershoots. Undershoots may appear if we include the nonlocal interaction 
into the model. That is the coupling between chains, e.g. the hydrodynamic interactions, 
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which are neglected in the present model. Another is that we suppose the shear flow is 
homogeneous, which in general is not. At last we can say that in polymer melts under shear 
non-uniform segmental orientation and stretch are hard to be observed directly in exper- 
iments, but their effect on rheological behavior can be detected, though the influences of 
the microscopic non-uniform behavior on them is not remarkable. However, the processes 
of crystallization and the dynamics of phase separation (Maurits and Fraaije, 1997) should 
be significantly affected by these contour dependent quantities, as the segmental dynamics 
plays a crucial rule in these processes. In fact the kinetic coefficients are related to the 
non-uniform features (Kawasaki and Sekimoto, 1988; Kawakatsu, 1997). It is of great sig- 
nificance to investigate the segmental dynamics and their affiliated properties theoretically 
by microscopic or mesoscopic models. 

ACKNOWLEDGMENTS 

We thank Dr. Bing Miao for the helpful discussion and suggestions. This work is sup- 
ported by National Natural Science Foundation of China (NSFC) 20490222, 20874111, 
50821062 and the Grant from Chinese Academy of Sciences KJCX2-YW-206. 



19 



APPENDIX A 



In this Appendix we describe the process and methods for the calculation of consti- 
tutive equations. The DE strain tensor without the independent alignment approxima- 
tion, i.e. Eq. (17), is integrated numerically. The unit vector of orientation is chosen as 
u = (sin 9 cos <p, sin 9 sin <p, cos 9) 



Q 



fE-uUE-u) 



■>■>) 



|E-u| / (|E-u|)o 

(U X + JUy) ■ Uy \ 



(Al) 



\ (1 + 2ju x u y + 7 2 ^)i/2 / Q ((1 + 2 1Ux u y + 7 2 «g) 1/2 )o ' 
where 7 is the shear deformation, (. . .) denotes the average over the isotropic state, i.e 
(. . .) = / du/47r(. . .). Then we have 

a{ 7 ) = {{l + 2 1 u x u y + 1 2 u 2 y ) l / 2 ) 

i-K r2ir 

— / sin^rf^ / dip(l + 27 sin 2 #sin if cos if + 7 2 sin 2 9 sin 2 y?) 1 ' 2 
47T Jo J 
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dx / dy 
Jo 



1 + 27 2 ( 1 ~ x2 ) + 7( 1 ~ x2 ) sin(47r?/) - ~ x2 ) cos(47ry) 



-,1/2 



(A2) 



Here, the variable transformations, i.e. x = cos 9, y = ip/2ir, are introduced. Then we have 

(1 — x 2 ) sin(47n/) + 7(1 — x 2 ) — 7(1 — x 2 ) cos(Any) 



Qx 



- / dx I dy 



a (7) 2 Jo Jo 



1 + |7 2 (1 — x 2 ) + 7(1 — x 2 ) sin(47r?/) — ^7 2 (1 — x 2 ) cos(47ry) 



i/2- 



(A3) 



With the same procedure, we obtain 

1 ,1 ,1 

Qxx Qyy / T / dx / dy 

a 7 Jo Jo 



x 



1 - I7 2 ) (1 - x 2 ) cos(47ry) + 7(1 - x 2 ) sin(47r?/) + |-y 2 (l - x 



-2- ( A4 ) 



1 + ^7 2 (1 — x 2 ) + 7(1 — x 2 ) sin(47ry) — |7 2 (1 — x 2 ) cos(47n/) 
The two-dimensional integrals of Eqs. (A2)-(A4) are integrated out using the extended 
Simpson's rule in each dimension, respectively. 

The constitutive equations under start-up of steady shear are solved numerically. We 
discretize Eq. (18) using the Crank-Nicolson scheme, and define F as 



F = u(i,j + 1) - q(i,j + 1) + u(i,j) + q(i,j), 



(A5) 
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where i and j denote the discretized contour variable grid points and time variable grid 
points, respectively. Here, u(i,j) is given by 

3Z(3 At 3 + a 4 q 4 (i,j) q(i,j) 
7r^ (AsJ^ [1 — a 2 q A {i,j)\ z 1 + q[i,j) 

3Z(3a At 4a 3 q 3 (t,j) + 12aq q(i,j) [q(i + 1, j) - q{i - 1, j)] 2 
vr 2 (As) 2 [l-aY^j)] 3 1 + <?(m) 4 

+-At 7 ^(z, J ) i + ^ (A6) 

where As and At are the contour and time grid steps, respectively. The contour grid step 
is a constant in our calculation, while there are two kinds of time step if the shear rate is 
higher than 10/7^. The grids density around the maximum value of Q X y{l) is larger than 
that the place far from it, so that the effects arising from the sharp changes of the value of 
Qxy{l) around the maximum value can be observed. There are N s — 1 coupled equations if 
we set F = for each time grid, and N s is the grid number of of the contour coordinate. 
The Newton-Raphson method is used to find their roots at each time grid. The number of 
time step is N t . In this paper, we choose As = 5 x and At has a range from 5 x 10~ 6 
to 4 x 10 -4 depending on the shear rate. For a larger shear rate, we choose a smaller At. 
Nt is determined by the time grid steps. Eq. (19) can be written as 



dG(s Q ,t) d 2 G{s ,t) dG(s ,t) 

= A (s , t) —y Ai(s , t) — h A 3 (s , t)G{s , t), 



dt 
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ds 
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where 
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A (s ,t) = l/(7rV), 

(do 1 ds 
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Eq. (A7) can be discretized as 
G(i,j + l)-G(iJ) 
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where 



( 



d 2 G 



dsl 




G(i+l,j)-2G(i,j) + G(i-l,j) 
(As) 2 



) 



G(i+l,j)-G(i-l,j) 



2As 



The complete procedure is as follows. For every time step, we first assume that S xy (i,j + 
1) = S xy (i,j). We solve the nonlinear equations of F = (F is given by Eq. (A5)), then 
we obtain q(i,j + 1). Insert these into Eq. (A8) we obtain G(i,j + 1). From Eq. (16) we 
obtain the new values of S xy (i,j + 1). By inserting these values back to Eq. (A5) and let 
F = 0, we obtain the new values of q(i,j + 1), then new G(i,j + 1), and subsequently 
Sxy{i,j + 1)- Repeat the above steps until the expected convergent value (10~ 12 ) is reached. 
The calculation of constitutive equations under steady shear is much easier than that under 
the start-up of steady shear, as in the former case the orientational tensor Q and strain q 
are time independent, and the calculation procedure are the same as in the case of start-up 
of steady shear at a certain time grid. 
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Figure Captions 

Fig. 1 Non-uniform segmental stretch with different shear rates under steady shear. 
7=0.1, f, fO, 100, 500, respectively, and Z = 20, a = 0. 

Fig. 2 Non-uniform segmental orientation with different shear rates under steady shear. 
7=0.001, 0.1, 1, 10, 100, 500, respectively, and Z = 20, a = 0. 

Fig. 3 Comparison of non-uniform segmental orientation of DEMG model and the present 
model under steady shear for a = and Z=20, 50, respectively. 

Fig. 4 Steady shear viscosity and shear stress for a = and Z=10, 20, 30, 50, respectively. 

Fig. 5 First normal stress differences and primary normal stress coefficients under steady 
shear for a = and Z=10, 20, 30, 50, respectively . 

Fig. 6 Comparison of the stresses and the first normal stress differences under steady 
shear between the present model (lines) and the experimental data (closed symbols). The 
parameters are chosen as Z = 20 and a = 0, 0.2, 0.4, 0.6, respectively. 

Fig. 7 Comparison of the stresses and the first normal stress differences under steady 
shear between the present model (lines) and the experimental data (closed symbols). The 
parameters are chosen as a = and Z = 10, 20, 30, 50, respectively. 

Fig. 8 The transient behavior of segmental stretch at Sq = for 7=0.1, 1, 10, 50, 100, 
respectively, and Z = 20, a = 0. 

Fig. 9 The transient behavior of segmental orientation at so = for 7=0.01, 0.1, 1, 10, 
100, 200, respectively, and Z = 20, a = 0. 

Fig. 10 The transient behavior of shear stress under different shear rates for 7=0.001, 
0.01, 0.1, 1, 10, 100, respectively, and Z = 20, a = 0. 

Fig. 11 The transient behavior of the first normal stress differences under different shear 
rates for 7=0.1, 1, 10, 100, 200, respectively, and Z = 20, a = 0. 

Fig. 12 Comparison of the shear viscosities with the experimental data (closed symbols). 
The parameters are chosen as Z = 20 and a = 0, 0.3, 0.6, respectively. 

Fig. 13 Comparison of the first normal stress differences with the experimental data 
(closed symbols). The parameters are chosen as the same as those in Fig. 12. 
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